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Abstract. We show that the shape of the observed distribution of Milky Way (MW) satellites is inconsistent with being drawn 
from a cosmological sub-structure population with a confidence of 99.5 per cent. Most of the MW satellites therefore cannot 
be related to dark-matter dominated satellites. 
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1. Introduction 

Calculations of structure formation within the framework of 
cold dark matter (CDM) cosmology show that Milky Way 
(MW) type systems have the same scaled theoretical dis- 
tribution of sub-haloes as rich galaxy clusters, and within 
500 kpc they should contain about 500 sub-haloes with masses 
M > 10 8 M (Moore et al. [T999l Klypin et al. [T999l Governato 
et al. 2004). However, only 13 dwarves have been found within 
a distance of 500 kpc around the MW. The observed dwarves 
may only sample a sub-set of the actually present CDM sub- 
structures (Stoehr et al. 2002; Hayashi et al. 2003; Bullock et 
al. 2000; Susa & Umemura 2004; Kravtsov, Gnedin & Klypin 
2004). Such biasing could be the result of complex early bary- 
onic physics that cannot, at present, be treated theoretically in 
sufficient detail, but Kazantzidis et al. (2004) point out that this 
cannot be the entire solution. 

An additional path to testing predictions of CDM cosmol- 
ogy is to compare the shape of the observed satellite distri- 
bution to the theoretical shapes (Zaritsky & Gonzalez 1999; 
Hartwick 2000; Sales & Lambas 2004). The sub-structures fall 
inwards from filaments that are thicker than the virialised re- 
gions of the hosts. However, within its virialised region, the 
number distribution of sub-structure in a theoretical host halo 
follows that of its dark-matter (DM) distribution. CDM models 
predict the host DM haloes to be oblate with flattening increas- 
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ing with increasing mass and radius (Combes 2002; Merrifield 
2002). The ratio of minor to major axis of the DM density dis- 
tribution has the value q& = 0.7 + 0.17 for MW-sized haloes 
within the virial radius. The intermediate to major axis ratio is 
q' d ^ 0.7 (Bullock 2002). When dissipational baryonic physics 
is taken into account the haloes become more axis-symmetric 
(larger q' d ) and more flattened, qa = 0.5 + 0.15 within the virial 
radius. The minor axis is co-linear to the angular momentum 
of the baryonic disk (Dubinski 1994). Prolate haloes do not 
emerge. The empirical evidence is that the MW dark halo is 
somewhat flattened (oblate) with q^ > 0.8 within R < 60 kpc 
(Oiling & Merrifield 2000, 2001; Ibata et al. 2001; Majewski et 
al. 2003; Martmez-Delgado et al. 2004). Beyond this distance 
the shape is likely to be more oblate (Bullock 2002), but invok- 
ing continuity the axis ratio q^ cannot change drastically. The 
theoretical sub-structure distribution of MW-type hosts must 
therefore be essentially isotropic (Ghigna et al. 1998; Zentner 
& Bullock l2003l Diemand, Moore & Stadel 2004; Kravtsov et 
al. 2004; Aubert et al. l2"U04l 

If the MW dwarves do indeed constitute the shining frac- 
tion of DM sub-structures, then their number-density distribu- 
tion should be consistent with an isotropic (i.e. spherical) or 
oblate power-law radial parent distribution. This is assumed to 
be the case by most researchers, given the relatively small num- 
ber of satellites. With this contribution we show that, despite 
its smallness, the MW satellite sample is inconsistent with a 
cosmological sub-structure population. We do this by concen- 
trating on the most elementary facts, namely purely on the po- 
sitions of the satellites. 
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Table 1. Dwarf galaxies within the vicinity of the MW. The first column is a running number used throughout this text; the 
parentheses contain the running dwarf number used in § after excluding SMC and UMi. D and eD are the distance and its 
uncertainty, respectively. /, b are the galactic longitude and latitude, respectively, as seen from the Sun and defined such that 
I — 0, b — points towards the Galactic centre which is assumed to lie at a distance D Q = 8.5 kpc from the Sun, and / increases 
in anticlockwise direction. The Galactocentric distance of the dwarf is given by R. The name of the dwarf is contained in the 7th 
column. The data are from Mateo ( 1998 table 2), except that for the LMC D and eD are taken from Salaris et al. (2003} and 
Clementini et al. (2003 1, and likewise for the SMC from Dolphin et al. |2001 1. The remaining columns contain the plane-fitting 
results for the innermost N dwarfs (§0: R cut is the largest distance to the Galactic centre of this sample, and the fitted plane has a 
root-mean square height A and a distance to the Galactic centre Dp. For comparison, the final column lists the root-mean square 
height A2(/?cut) for samples of 4 x 10 5 theoretical dwarves with an isotropic isothermal radial number density profile (p = 2) and 
radial cutoff R cut . 
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2. Dwarf galaxies near the Milky Way 



Table fTJlists distances and coordinates of the N = 16 dwarves 
closest to the MW. Given these data Galactocentric coordi- 
nates are calculated, Xp, = -D Q + D sin(90° - b) cos(Z), Yp, = 
Dsin(90°-fr)sin(0,Z D = Dcos(90°-fr), with uncertainties de- 
rived from the uncertainties in D. The data are compared with 
an isotropic power-law density distribution, p(R) = p R p , 

where R = -yJx D + + Z Q . Using the Kolmogorov-Smirnov 
(KS) test it can be shown that the cumulative dwarf sample 
is consistent with a radial near-isothermal density distribution, 
the 5 per cent confidence margin being 1.8 < p < 2.6 for the 
N = 1 1 innermost dwarves. The solutions shift to larger p as 
further outlying dwarves are added. This is a similar behaviour 
as seen in theoretical distributions of sub-haloes (e.g. fig. 5 in 
Zentner & Bullock l2003l . 

The data are plotted in Fig. [Rafter clockwise rotation by an 
angle <p = 167°. 9 about the Z-axis, X = X D cos0 + Y D sin</>, Y = 
-Xp, sin<p + Yd cos</> and likewise for the uncertainties. The dis- 
tribution is highly anisotropic and planar in form. It is the 
aim of this contribution to quantify the significance of this 
anisotropy. A rotation of Fig. ^ by 90° shows the distribution 
to be approximately disk-like (Fig. 13 . 



3. The satellite plane 

A plane can be described by the Hesse form, n • (x — p) = 0, 
where n is the normal vector, p a vector pointing from the ori- 
gin (the centre of the MW) to a point in the plane, and x an ar- 
bitrary vector from the origin to the plane. With n = (n\ , ri2, «3) 
and x(i) = (X(i), Y(i),Z(i)) being the coordinates of the galax- 
ies, d(i) = n\X{i) + rt2Y(i) + n^,Z(i) - Dp, becomes iden- 
tical to Hesse form if d(i) = 0; d(i) being the distance of 
the z'th dwarf to the plane. Dp = n • p is the shortest dis- 
tance of the plane to the origin. The problem of finding the 
plane can thus be reduced to a least-squares linear regres- 
sion problem, where the aim is to find the coefficients, n,, Dp 
with the condition J^ =i nf = 1, that minimise £]j=i d 2 (i). To 
achieve this the method of normal equations using Gauss- 
Jordan elimination is employed to solve the set of linear equa- 
tions (Press et al. 1992 1. For each fitted plane the root-mean 
square height of the resulting disk distribution is calculated, 

ACR cut ) = -yJ(l/N) d 2 (i). Note that the applied minimisa- 
tion does not include the location of the Galactic centre as a 
constraint. Thus, in principle the fitted plane to a small number 
of dwarves (N < 12) could lie far away from the Galactic cen- 
tre. The weights that do enter the regression are merely given 
by the uncertainties in distance. The direction of the normal 
vector, or the location of the pole of the plane, lp, bp, follows 
from = arccos(«3), b' p = 90° - 6,l' v = arcsin((«2)/sin(6')). As 
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Fig. 1. The position of the innermost 1 1 MW satellites (Tabled 
as viewed from a point located at infinity and I = 167°. 91. The 
MW disk is indicated by the horizontal line -25 < X/pc < 
25, and the centre of the coordinate system lies at the Galactic 
centre. The dashed line marks the fitted plane for 7Y = 1 1 seen 
edge-on in this projection. 
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Fig. 2. As fig.^but viewed from / = 77°.91. The fitted plane is 
here seen face-on. 
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Fig. 3. The position on the Galactic sky of the poles of the 
planes fitted to the dwarves of Table [l] Plotted are bp = —b' p 
and lp — /p + 180° and the number of dwarves used for the fit 
ranges from N = 16 down to N = 3 (TableQ. The cases for 
N = 3,4,5 are indicated with numbers. The others cluster very 
tightly around l P a 168°, bp ~ -16°. The likely position of the 
orbital poles of the LMC, SMC, Draco and UMi are indicated 
by the solid curves (from fig. 3 in Palma et al. 2002). 



no kinematical information is included the direction of the pole 
is ambiguous, bp = b' p , lp = l' p , or bp = —b' p , lp - l' p + 180°. 

Tabled lists some results of the fitted plane for a decreas- 
ing number of dwarves. The empirical disk height, A, is always 
much smaller than the theoretical height, A2, for an isothermal 
and isotropic model number density distribution centred on the 
origin of the MW. The MW dwarves thus appear to be dis- 
tributed as a great disk with a ratio of height to radius < 0.15. 

The poles of the planes and the orbital poles of the dwarves 
LMC, SMC, Draco and UMi agree remarkably well (Fig. |3}. 
This is surprising because the results are obtained using com- 
pletely different methods. The position of the poles of the 
planes found here depend only on the spatial distribution of 
the dwarves. In contrast, an orbital pole is the direction of the 
orbital angular momentum and relies on the direction of the 
measured proper motion of the respective object. Sgr is on a 
polar orbit but has a kinematical pole (I « 280, b » 0, Palma, 
Majewski & Johnston 2002 1 lying approximately at a right an- 
gle to the great disk and to the MW disk. On the basis of the 
weakly bound core of Sgr which makes it difficult for Sgr to 
survive the many orbits implied by its current angular momen- 



tum, Zhao (1998) proposed that it may have been scattered 
into its present low-pericenter orbit by an encounter with the 
LMC/SMC about 2-3 Gyr ago. Sgr contributes the most de- 
viant cos(w) value in the sample because it is closest to the 
MW centre and thus high above the local great disk. Taking 
Sgr out of the sample would increase the discrepancy, quanti- 
fied in § |3 between the dwarf sample and the hypothesis that 
they are the visible cosmological sub-halo population. 

4. The likelihood 

The null hypothesis is that the N observed dwarves are drawn 
from a cosmological population. We therefore need to establish 
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the probability that the observed distribution is drawn from a 
spherical parent distribution. 

The vector pointing from the Galactic centre to the closest 
point, P c i, on the plane is dp = D P ra, and the vector from this 
point P c i to a dwarf is x' = -dp + x. The angle, a>, between the 
normal vector and the dwarf as viewed from P c i is then given 
by cos(w) = n • x'/\x'\. The cumulative distribution of cos(<y) 
about the fitted plane is calculated for the observed sample us- 
ing the innermost N dwarves, and also for N m - 10 5 model 
dwarves distributed according to the theoretical parent radial 
power-law distribution which is centred on and isotropic about 
the Galactic centre. The KS test quantifies the confidence that 
can be placed in the hypothesis that the observed sample stems 
from this parent distribution. The results, plotted in Fig.0] show 
that this hypothesis can be rejected with a confidence of better 
than 98 per cent, and even 99.6 per cent for N > 11. This comes 
about because the real sample is deficient near the poles of the 
great disk. 

Orbital pole analyses have been showing that the SMC, 
UMi and the LMC form a kinematical family (Palma et 
al. 2002). Taking these two objects out of the sample, 
kinematically-linked dwarves are removed with the expectation 
that the remaining dwarves should be more consistent with the 
isotropic parent distribution. As the thin curves in Fig.0]show 
this is not the case. Instead, the probabilities that the N = (9) 
sample without the SMC and UMi stems from an isotropic par- 
ent distribution is reduced (as compared to the N - 9 sam- 
ple). This comes about because the two dwarves are relatively 
close to the Galactic centre thus adding relatively large u> an- 
gles when they are included. 

The disk-like distribution of the dwarves lying nearby to 
the MW noted in Fig.^is therefore highly significant. The lo- 
cal dwarves do not stem from an isotropic distribution. Their 
distribution is therefore severely at odds with the sphericity of 
the MW dark matter halo, and even more at odds with an oblate 
halo having the same orientation as the MW disk. 

5. Concluding remarks 

Cosmological models can be tested among other ways by 
comparing the theoretical sub-structure distribution with ob- 
served satellite distributions. The theoretical distribution con- 
tains about 500 sub-haloes within approximately 500 kpc of 
a MW-type galaxy and follows an approximately power-law 
radial distribution with p < 2 and is essentially isotropic. The 
well-known MW distribution contains only a dozen dwarves, 
is indeed consistent with the theoretical radial distribution but 
is highly anisotropic. The anisotropy is such that the MW 
dwarves form a disk-like structure with a root-mean-square 
height of 10 - 30 kpc which lies nearly perpendicularly to the 
plane of the MW. The pole of this great disk lies close to the or- 
bital poles of the LMC, the SMC, Draco and Ursa Minor. The 
distance of closest approach of the plane to the Galactic cen- 
tre, Dp ^ 2 kpc, is much smaller than the radial extent of the 
Galactic disk 20 kpc) or even the root-mean square height, 
A, of the disk of satellites (Dp <K A). This is a strong indication 
that the sample of dwarves within about 250 kpc is relaxed in 
the Galactic potential. Their orbits must be confined within the 




P 



Fig. 4. The probability, P C osw, that the observed dwarf sample 
stems from a parent isotropic radial power-law density distri- 
bution with index p. The number of innermost dwarves in the 
sample (Table [3 is indicated by the numbers. Thus 11, for 
example, means that the innermost N = 1 1 dwarves (out to 
and including Leol) are compared with the isotropic power- 
law distribution. The thin curves (and numbers in parentheses) 
are probabilities calculated by excluding SMC and UMi from 
the data; here N = (9) incorporates all dwarves except SMC 
and UMi out to and including Leol (Table^. The increase of 
P C0S0J with decreasing N and (AO is a result of weakening confi- 
dence as the number of observed data is reduced. The decrease 
of Pcosw for p > 2.4 comes from the theoretical distribution be- 
ing increasingly concentrated towards the Galactic centre while 
the plane lies off-centre (Dp > 0). 

great disk because the likelihood of obtaining such a disk-like 
dwarf distribution given a true underlying isotropic distribu- 
tion (that ought to match the sphericity of the MW DM halo) is 
less than 0.5 per cent. This result persists even after removing 
the kinematically related SMC and UMi from the analysis. A 
distribution of polar orbits with arbitrary eccentricities and ori- 
entation of orbital planes is also excluded with the same confi- 
dence because it leads to an isotropic distribution of dwarves. 
An oblate MW dark matter halo would yield an even larger 
discrepancy with the disk of satellites. 

An alternative approach is taken by Hartwick (2000) who 
argues that the 10 satellites within 400 kpc (the LMC and SMC 
are combined to one satellite) map the MW DM halo shape 
and form a highly inclined and highly prolate system with mi- 
nor/major axis ratio q& « 0.03 - 0.05. However, the extreme 
triaxiality derived in this way is completely inconsistent with 
the observational and theoretical shapes of CDM host-haloes 
and sub-structure distributions (§[D- 
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The approach taken here differs by noting the very sig- 
nificant mismatch between (i) the disk-like satellite distribu- 
tion, (ii) the independent empirical constraints on the shape of 
the MW dark matter halo, and (iii) the theoretical shapes of 
CDM host haloes (§[D. In the view presented here, the mis- 
match between the number and spatial distribution of MW 
dwarves compared to the theoretical distribution challenges the 
claim that the MW dwarves are cosmological sub-structures 
that ought to populate the MW halo. 

A more natural and more conservative (by not resorting 
to exotic physics) explanation for the MW dwarf distribution 
in a great disk with a ratio of height to radius of 0.1-0.2 
would appear to be in terms of a causal connection between 
most of them. This could be the case if most of the dwarves 
stem from one initial gas-rich parent satellite on an eccen- 
tric near-polar orbit that interacted with the young MW, per- 
haps a number of times, forming tidal arms semi-periodically 
as its orbit shrank. The early gas-rich tidal arms may have 
condensed in regions to tidal dwarf galaxies, as is observed 
in present-day interacting gas-rich galaxies (e.g. Knierman et 
al. 2003; Weilbacher, Due & Fritze-v. Alvensleben 2003). The 
LMC may be the most massive remnant of this larger satel- 
lite, while the lesser dwarves may be its old children (Lynden- 
Bell 1976 1. The Magellanic Stream may be just such a newly 
formed but meagre tidal feature (Kunkel 1979 1, and the align- 
ment of the disk of satellites with the surrounding matter distri- 
bution (Hartwick 2000) may simply result from the gas-rich 
parent satellite coming-in from that direction. The different 
chemical enrichment and star-formation histories of the vari- 
ous dwarves (e.g. Ikuta & Arimoto 2002, Grebel et al. 2003 1 
may in this case be a result of their different initial masses that 
will have been significantly larger than their present-day bary- 
onic masses (Kroupa 1997 1 and the complex interplay between 
stellar evolution, tides, gaseous stripping and gas accretion dur- 
ing the orbits within the MW halo, none of which are presently 
understood in much detail. The simulations of Kroupa have 
shown that ancient tidal-dwarf galaxies may appear similar to 
some of the observed dSph satellites. 

The sub-structure under-abundance problem extends to fos- 
sil galaxy groups where early photo-evaporation could not have 
removed baryons from the sub-structures (D'Onghia & Lake 
2004), and a sub-structure-overabundance is evident for rich 
clusters (Diemand et al. 2004). CDM cosmology thus faces a 
sub-structure challenge on all mass scales. 
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